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We analyze the photon statistics of a weakly driven optomechanical system and discuss the ef- 
fect of photon blockade under single photon strong coupling conditions. We present an intuitive 
interpretation of this effect in terms of displaced oscillator states and derive analytic expressions for 
the cavity excitation spectrum and the two photon correlation function g^^\0). Our results predict 
the appearance of non-classical photon correlations in the combined strong coupling and sideband 
resolved regime, and provide a first detailed understanding of photon-photon interactions in strong 
coupling optomechanics. 
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The implementation of strong optical non-linearities on 
a single photon level is one of the central goals in quan- 
tum optics with a significant practical relevance for ap- 
plications ranging from optical computation [Tj to quan- 
tum information processing [2] and photonic quantum 
simulation schemes [3]. The prototype system that has 
been widely studied in this context is cavity QED [4^ 
where under strong coupling conditions effective photon 
non-linearities result from the hybridization between the 
optical field and a single atom. Recently, a fundamen- 
tally different type of light-matter interaction has at- 
tracted a lot of attention, which is the radiation pressure 
coupling between light and mechanical motion studied 
in optomechanical systems (OMS) [S]. In most exper- 
iments today radiation pressure forces are fairly weak 
and non-linear optical effects [7] occur in the classi- 
cal, high photon number regime, where enhanced but 
linear photon-phonon interactions [SI |H] are investigated 
for cooling [10] or the preparation of mechanical quantum 
states fll'. However, strong optomechanical interactions 
with single photons in analogy to cavity QED, are within 
reach of new generations of nano-fabricated OMS [T^] or 
superconducting devices [S] and are already nowadays 
accessible in analogous cold atom experiments [13l [M] . 
This could open up a new route towards non-linear quan- 
tum optics, which avoids single atom strong coupling and 
where instead photon-photon interactions are mediated 
by the motion of a macroscopic object. 

In this work we study OMS in the regime where the 
single photon coupling go is comparable with the cavity 
decay rate k. In contrast to previous studies [TSHTO] we 
here focus explicitly on the consequences of strong cou- 
pling for the quantum statistics of light, with the aim to 
identify the mechanism for photon-photon interactions in 
this system and the experimental conditions under which 
such effects could be observed in experiments. To do so 
we consider a weakly driven OMS as shown in Fig. [l] 
and evaluate the two photon correlation function g^'^\Q). 
This quantity provides a direct experimental measure 
for non-classical anti-bunching effects, i.e. (7^^^(0) < 1, 
and the limit ^'•^-'(O) — ^ indicates a complete photon 
blockade pOH22] . where strong photon-photon interac- 
tions prevent multiple photons from passing through the 




FIG. 1. Setup for the detection of photon blockade effects in 
optomechanical systems (OMS). The OMS is weakly excited 
by a coherent laser field and the statistics of the output field 
is inferred from photon coincidence measurements. 



cavity at the same time. Our analysis shows that apart 
from 5o signatures of non-classical light in OMS depend 
crucially on the relation between the cavity linewidth n 
and the mechanical frequency LOm and appear only un- 
der quite stringent conditions k < gQ,uJm- However, this 
regime is within reach of experiments ll2H14j where 
the observation of photon blockade would provide the es- 
sential ingredient for potential applications of OMS as a 
quantum non-linear device. 

Model. We consider a setup as shown in Fig. [T] where 
the frequency of an optical cavity mode is modulated by 
the linear displacement of a mechanical oscillator. The 
cavity is excited by a weak coherent laser field and the 
photon statistics of the transmitted light is analyzed us- 
ing standard photon counting techniques |22j . In a frame 
rotating with the laser frequency uj^ the Hamiltonian for 
the OMS is {h = 1) 



Hop = i^,„-Aoc^c^-^5fc(4 + &fc)c^c+^£ {c'' - c) , (1) 
fc 



where c is the bosonic operator for the cavity mode, £ 
is the driving strength and Ap = uj^ — lOq is the de- 
tuning of the laser from the bare cavity frequency lu^. 
The bosonic operators bk represent the mechanical eigen- 
modes of system which evolve under the free Hamilto- 
nian H„i = X^fe'^fe^I^fc a^'i couple to the cavity with 
a strength g^- In a general device the bk's account for 
different vibrational modes of the resonator as well as 
mechanical modes of the support and the system spe- 
cific details of the OM interactions are summarized by 
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FIG. 2. Level diagram of the isolated, single mode OMS where 
T] — go/c^m and Ag = (jg/i^m- For a photon number state 
\nc) the radiation pressure displaces the resonator equilibrium 
ricXrj. As a result the energy of the photon states is lowered 
by X Ag and leads to different resonance conditions for the 
first and the second laser photon exciting the cavity. 



the spectral density J{uj) = ^J2k9k^(^~^k) For 
concreteness we will below focus explicitly on the case 

which models a mechanical mode of frequency cj„i cou- 
pled to an Ohmic bath. The dimensionless parameter 
r] = go /ujm is chosen such that in the limit of a high me- 
chanical quality factor Q we recover the standard model 
for a single mode OMS f5] with a coupling constant go. 

The cavity field is coupled to the electromagnetic vac- 
uum modes of the environment and in the limit of a single 
sided cavity we model the resulting dissipative dynamics 
by a quantum Langevin equation 

C(t) = t[Hop, c{t)] - Kc{t) - V2KU,{t). (3) 

Here k is the cavity field decay rate and fin{t) a ^- 
correlated noise operator. Photon counting and pho- 
ton coincidence measurements of the cavity output field 
fout{t) — fin{t) + V'2'Kc{t) provide information about the 
photon statistics of the cavity field [22] . 

Displaced oscillator states. For the observation of pho- 
ton blockade effects we are interested in the regime of low 
photon numbers where the driving field £ only weakly 
perturbs the OMS. Therefore, to proceed it is conve- 
nient to change to a displaced oscillator representation. 
Hop — >■ UHopU\ which diagonalizes Hop in the limit 
£ — >■ and is defined by the unitary transformation 
U = e-'P'"^" and P ^ iJ2k(9k/^k){bl - hk). We obtain 

Hop = Hrn - Ac^c-AgC^c^cc + i£ {c^e-'^ ~ e'^c) , (4) 

where we have introduced a shifted detuning A = Ag + 
Ag and a photon non-linearity Ag = f /q dujJ{Lo)/uj, 
where An — g'^/ujm for the single mode model defined 
in Eq. pi). The origin of this effective photon-photon 
interaction can be understood from the fact that in an 
isolated system the radiation pressure force displaces the 
resonator equilibrium by an amount proportional to the 



photon number ric and thereby lowers the energy of this 
photon state by x Ag. This is illustrated in more detail 
in Fig. [2] and already explains the basic mechanism for 
photon blockade. If the driving laser is on resonance 
with the |0c) — > |lc) transition, i.e. A = 0, the same 
|lc) — > |2c) transition is detuned by 2Ag and will be 
suppressed for Ag > n. However, this simple picture is 
based on the level structure of the isolated OMS [13 US] 
only and ignores phonon sideband transitions and other 
dynamical aspects of the problem which will be addressed 
by the following more rigorous analysis. 

Excitation spectrum. We first study the cavity excita- 
tion spectrum iS'(Ao) := \Ymt^oo{c\t)c{t)) /tiq, which is 
normalized to the resonant photon number of the unper- 
turbed cavity, = /k? . The OMS is initially prepared 
in the state p(0) = |0c)(0c| (8) Pth where pth is the ther- 
mal equilibrium state of the mechanical modes. For a 
finite, but weak driving field the dominant contribution 
for S'(Ao) arises from terms in the Heisenberg operator 
c{t) which are linear in £ and from the displaced oscilla- 
tor representation of Eq. ([s]) we obtain 

c{t) = {iA~K)c{t) + e-'P^'\£~V2Kf.Ut)) + 0{£^). (5) 

On the same level of accuracy the operator P{t) 
can be approximated by the free evolution P{t) = 
^-iHr,^tp^iHmt after integrating Eq. ([5| we find 

S{An)^KRe dr e'^'^-^^^ e'^^^'^ + 0{£). (6) 
Jq 

Here e"^^^'^-' — (e'^'-'^'e"*^^*'-') is the equilibrium cor- 
relation function of the mechanical displacement opera- 
tor For later convenience we define /2 = /2(w,r) = 
2„g-ii.jT g^j^j write the function F2{t) in the general form 

mn}) = - Hdu:^ [(7VM+i)/,+7VH/,!] , 

TT Jo ^ 

(7) 

where N{uj) = l/(e''"/'^^-^ — 1) is the bosonic equilibrium 
occupation number for a support temperature T. 

Fig. [3] shows the cavity excitation spectrum S'(Ao) 
for two different values of K/cOm and for J(w) defined 
in Eq. In the bad cavity limit k ^ uJm we can 

approximate F2{t) ~ i^g''' + '''^/(4T'^) and we iden- 
tify an optomechanical dephasing mechanism with a 
timescale T'^ = 2go^/N + 1/2 where N := iV(a;„). 
This leads to a broadening of the spectrum and a grad- 
ual change from a Lorentzian to a Gaussian lineshape 
5(Ao) ~ y^nT^^e"^"'^^ for very large values of go. A 
completely different behavior is found in the sideband 
resolved regime k <^ ujrn- Here we observe a red-shift 
of the zero phonon line (ZPL) towards Aq = — Ag and 
the appearance of additional resonances at multiples of 
the mechanical frequency uJm- These peaks result from 
phonon assisted excitation processes. 

For a more detailed quantitative discussion of S'(Ao) 
we now focus on the limit Q ^ 1 of a weakly damped 
mechanical mode. In this regime F2{t) = -FJ (''') + ^-^2 ('''): 
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FIG. 3. Cavity excitation spectrum S'(Ao) for different values 
of the coupling parameter rj = go/u)m and a) Hi/ujm = 4 and 
b) K/um = 0.1. In b) the zero phonon line (ZPL) indicates 
the position of the phonon number conserving transition. In 
both plots T = and Q = 150. 



where — rf svi\{uj„iT)e~i^ and F2{t) ~ rT + 77^(2A^ + 
1)(1— cos(a;„ir)e~ 2"^). Here 7 = ijJmlQ is the mechanical 
damping rate and V is an additional decoherence rate 
which arises from the low frequency part of J(w). It 
vanishes for T — ^ and is given by F ~ 77^(2A^ + 1)7 for 
temperatures T > hLUm/kg. The approximate analytic 
expression for F2{t) allows us to expand the two point 
correlation function in Eq. ([6| and evaluate the integral 
over T [Ml. We obtain 



+ (Ao + Ag-nw„)2' 



(8) 



where A„ == e^" (2A'+i)/„(2r;V^(^+l)) {{N+l)/Ny 
and In{x) is the n-th order modified Bessel function. 
This result is familiar from the standard Huang-Rhys 
theory of phonon assisted excitation processes P5 ^'^d 
the positions and weights of the resonances can be un- 
derstood from different multi-phonon sidebands of the 
|0c) — >■ |lc) transition shown in Fig. [2] Apart from pho- 
ton loss the resonances are broadened by the mechanical 
decoherence rates where k„ ~ k + F + |n|7/2 for TV < 1 
and K„ ~ K -f 2F in the hight temperature limit. We 
emphasize that at T = both the appearance of phonon 
sidebands for k < ujm as well as the broadening of the 
cavity resonance in the opposite regime k > uj,n are pure 
quantum effects and provide a clear indication for single 
photon strong coupling optomechanics. 

Photon correlations. For weak driving the excita- 
tion spectrum is dominated by single photon events 
and does not contain information about photon-photon 
interactions. To proceed we now concentrate on the 
normalized equal time correlation function g^^^O) '■= 
limt^ao {c'c'cc){t)/{c^c)'^{t), where in addition to S'(Ao) 
we must evaluate the two photon correlation G'-^' := 
limt_j.oo(c^^(i)c2(i))/riQ. Following the same arguments 
as above we obtain up to the lowest relevant order in £, 

c\t) = 2{iA + iAg-Ky{t)+£e-'^'^'^c{t)+0{£^), (9) 



where we have already omitted an irrelevant noise term 
~ fin(t)- Together with Eq. ([s]) we finally obtain 



G(2) = 2k^ Re / dTi 



(10) 



^g2(iA+iAp-K)rig(-iA-K)r2g(iA-K)r3g-F4(ri,T2,T3) 



where e"^''^^'^'^) 



jP(0)„-»P(-r3)\ 



a four point correlation function of the mechanical dis- 
placement operator. The function Fji{{Ti}) can be ex- 
pressed in terms of Eq. ([t]) by setting /4 = fi{u),{Ti\) = 



(1 



^(1 



Since a general discussion of Eq. ( 10 ) is quite in 



volved we will from now on concentrate on the most 
relevant regime where the mechanical decoherence rates 
F and 7 can be neglected compared to k. How- 
ever, we first point out that in the bad cavity limit 
we can approximate F4({Ti}) « iAp(4Ti — T2 + T3) -|- 
(2ti - T2 + T3)2 /(4r2). Then, for Aq < T'^ we ob- 
tain g'^^O) ~ e^-^o^-^^VlViTrKr^) > 1 and we conclude 
that even for strong coupling 50 the photon statistics of 
an OMS in the bad cavity limit remains classical. 

Let us now consider the limit Q — ^ 00 where -Fl^dri}) ~ 
rf fi{ijJrm{Ti\) and as above, we use a series expansion 
of the correlation function in Eq. ( 10 1 to evaluate the 
integrals over the r^. We obtain 

X ^ Br) _m._ri 



2^3 



(11) 



(2k - i{2A+2Ag-puj,n)) ' 

where the coefficients Bn^m,p depend on the coupling pa- 
rameter 77 and follow from, 

p-')^/4(t^m,{Ti}) _ r> iUm{T2n-T3m-Tip) 



n,m,p 



(12) 



For example, for T — Q explicit expressions are given 
by Bn^ni,p = e-'^'i'' {ri^YWn.,p{ri)W,n.p{ri) /n\m\p\ where 
Wn,p{ri) = (-l)"U[-n, 1 - n -I- p, 77^] and U[a, &, x] is a 
confluent hypergeometric function. 

Discussion. In Fig. |4] a) we use Eq. ^ and Eq. (11) 
to evaluate g^^-'(O) and plot the result as a function of 
Ao and for k <C a;^. We observe a sequence of bunch- 
ing and anti-bunching resonances which can be qualita- 
tively understood from the level diagram shown in Fig. [2] 
where depending on Aq either the one or the two photon 
transition becomes resonant with different phonon side- 
bands. For a better understanding of this process we now 
consider the regime 77 < 1 and assume that the laser is 
tuned close to the ZPL of the one photon transition, i.e. 
A = Ag -I- A„ <C Then, still assuming k < w^, the 



dominant contributions to the sum in Eq. ( 11 ) arise from 
the terms n — m = Q and p = 0, 1 and combined with 
the n = terms in Eq. ([8| we obtain 

Go(k2 + A2) ri'^Ciin' ^ a-') 



'■ + {A+Ag)^ K^ + {A+Ag-Lj,n/2y 



(13) 
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FIG. 4. a) Dependence of (;''^'(0) on the bare laser detuning 
Ao for go/ujrn ~ 0.5 and n/iUm ~ 0.15. The dashed hne 
indicates the approximate result given in Eq. ( |13[ ). b) The 
minimum of 3'^' (0) with respect to Aq is plotted for different 
values of k and go- In both plots T = and Q 00. 



Here Co = Bq.o.oMo ^^id C*! = Bq 
zero temperature Cq = Ci = 1. Fi 
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ig- El 



77 Aq) and for 
a) shows that 

Eq. ( 13 1 provides indeed an excellent approximation of 
the first anti-bunching tip around Aq — — Ag and allows 
us to make the following analytic predictions. First, for 
K > go the minimum for g'-^^\0) occurs at A w k and 
scales as min{(;(^^(0)} « 1 — / {lo„iI^) . Therefore, as ex- 
pected, no significant anti-bunching effects appear unless 
the strong coupling condition go > ^ is achieved. In this 
regime the minimum of (7^^^(0) occurs at A = and 



min{g(2)(0)} 



(1 - 2772 



This result demonstrates that OMS can indeed exhibit 
a strong photon blockade effect where for 77 <C 1 the 
suppression of two photon events scales with the param- 
eter K^ujf^/gQ. However, rather than improving mono- 
tonically with increasing coupling strength the blockade 
reaches a minimum at 77 ~ 0.5 with a value g^'^^O) ~ 
20 X (k/lli„j)2. This minimum is a consequence of the 
|lc) — )■ 1 2c) transition getting into resonance with the first 
phonon sideband which occurs for rj — l/\/2- Therefore, 
the fidelity of the photon blockade effect is ultimately 
limited by the sideband parameter n/iOm and the effect 
vanishes as this parameter approaches one. 

Finally, we use Eqs. (|8| and (11) to evaluate the min- 
imum of (7^2) (0) for a large range of parameters go and k 
numerically. The results are plotted in Fig.|4]b) and show 
a clear boundary at k ~ Wm which - quite independently 
of the value of go - separates the regimes of classical and 
quantum correlations. We also see that in the sideband 
resolved regime n < Um the photon blockade exhibits a 
repetitive pattern and for g^ > uiml"^ no significant fur- 
ther improvement is achieved. These results show that 
the approximate result ( 13 1 already captures the essence 
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of the two photon blockade effect in OMS. 

In summary we have identified the mechanism for 
strong photon-photon interactions in OMS and studied 
the dependence of the two photon blockade effect on the 
relevant parameters go-, k, and iOm- Our results provide a 
guideline for future experiments and a first detailed the- 
oretical description of the two photon physics which is 
relevant for applications of OMS in the context of quan- 
tum information processing or quantum simulation. 
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